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Abstract. We consider that no mean magnetic field ex- 
ists during this epoch, but that there is a mean mag- 
netic energy associated with large-scale magnetic inhomo- 
geneities. We study the evolution of these inhomogeneities 
and their influence on the large scale density structure, by 
introducing linear perturbations in Maxwell equations, the 
conservation of momentum-energy equation, and in Ein- 
stein field equations. The primordial magnetic field struc- 
ture is time independent in the linear approximation, only 
being diluted by the general expansion, so that BR^ is 
conserved in comoving coordinates. Magnetic fields have a 
strong influence on the formation of large-scale structure. 
Firstly, relatively low fields are able to generate density 
structures even if they were inexistent at earlier times. 
Second, magnetic fields act anisotropically more recently, 
modifying the evolution of individual density clouds. Mag- 
netic flux tubes have a tendency to concentrate photons 
in filamentary patterns. 



Key words: Magnetohydrodynamics (MHD) - Relativity 
" Cosmology: large-scale structure of universe 

1. Introduction 

Recent measurements of intergalactic magnetic fields (see 
the review by Kronberg, 1994, and references therein) have 
provided evidence of the following facts: 

1. Magnetic fields of the order of are not uncom- 
mon. Values larger than IfiG have been found in all 
reported measurements (see also Feretti et al. 1995). 
Fields of this strength are quantitatively important be- 
cause they correspond to an energy density equal to 
that of CMBR. Kronberg (1994) even suggests that 
S/iG fields are ubiquitous. In this "background 
magnetism" would be in equipartition of energy with 
the CMBR. This possibility, even whilst theoretically 



attractive, is still based on a limited number of mea- 
surements and therefore will not be assumed here. 
2. Large magnetic fields have also been found in proto- 
galactic clouds (e.g. Wolfe, Lanzetta, & Oren 1992; 
Welter, Perry, & Kronberg, 1984). Probably, primor- 
dial magnetic fields to a large degree contributed to 
the present intergalactic ones. Kronberg suggested that 
magnetic fields play a role in the formation of struc- 
tures in the Universe. 

These facts, even if we do not know exactly how ubiq- 
uitous and how persistent in time intergalactic magnetic 
fields are, stimulate the analysis of their evolution and 
interrelation with density inhomogeneities. 

In this paper we consider a universe dominated by rel- 
ativistic particles. More precisely we have in mind a uni- 
verse dominated by photons before Equality. The equa- 
tions however are valid for any kind of dominant relativis- 
tic particles including hot dark matter. 

The study of the evolution of density inhomogeneities 
for a universe with photons and barions, when no mag- 
netic fields are present, is a classical topic (Weinberg 1972; 
Peebles 1980; Kolb & Turner 1990; Borner 1988; Battaner 
1996). It may be divided into three periods: 

i) Post-Recombination era. During this era a New- 
tonian analysis is appropriate, but nonlinear effects 
require rather sophisticated numerical techniques. 
Inhomogeneities grow as R first, becoming propor- 
tional to R^ and R^ when nonlinear effects become 
more and more important. The inclusion of mag- 
netic fields in the study of this epoch has been car- 
ried out by Wasserman (1978) , Coles (1992) and 
Kim, Olinto, & Rosner (1994), this last work in- 
cluding nonlinear effects. Classical treatments of 
Birkeland currents in the plasma Universe have 
been carried out by Peratt (1988) 

ii) Acoustic era. A relativistic treatment is neces- 
sary; viscosity and heat conduction must be in- 
cluded (Field 1969; Weinberg 1972; and others) 



2 



E. Battaner et al.: Magnetic fields and large scale structure in a hot universe. I. General equations 



as these effects explain the Silk mass. Inhomo- 
geneities do not increase during this era, which ends 
at Recombination, its beginning being dependent 
on the rest mass of the primordial cloud, around 
R = 10~^i?o- As far as we know, no attempt at 
introducing magnetic fields into this analysis has 
been made. 

iii) Radiation dominated era. This era ends when 
the acoustic one begins, and is therefore not per- 
fectly defined, roughly at i? ~ 10~^i?o, so it corre- 
sponds to a photon dominated universe. The begin- 
ning of this era is also rather indeterminate. During 
this epoch non-magnetic inhomogeneities increase 
as 

The inclusion of magnetic fields in the study of this third 
era is the objective of this paper. Basically, our objective in 
this paper is to extend the work by Wasserman (1978) and 
Kim, Olinto and Rosner (1994) to the radiation dominated 
era. We will deal with the evolution of magnetic fields and 
their influence on density inhomogeneities in a radiation 
dominated Universe. The mathematical procedure must 
be relativistic, but the inclusion of nonlinear and imper- 
fect fluid effects is not necessary, which greatly simplifies 
the problem. The upper time boundary will be placed at 
approximately 10~^i?0) before the Acoustic epoch, and 
near Equality. The lower time boundary is undefined, but 
in particular we consider a Post- Annihilation era, in order 
to avoid sudden jumps in the temperature of photons and 
because positron-electron and quark-gluon plasmas, which 
constitute the plasma state at earlier epochs, require an- 
other analytical formulation. Therefore, the period under 
study broadly extends from Annihilation to Equality. 

We consider that the evolution of magnetic fields is 
not perturbed by creation and loss processes. Some mech- 
anisms have been invoked for later stages ( Rees 1987; 
Lesch & Chiba 1995; Ruzmaikin, Sokoloff, & Shukurov 
1989; and others) but these probably do not affect the 
epoch studied here. Some mechanisms producing primor- 
dial fields, prior to Annihilation, are implicitly assumed 
(see Turner & Widrow 1988; Quashnock, Loeb, & Spergel 
1989; Vachaspati 1991; Ratra 1992; Enqvist & Olesen 
1993, 1994; Davis & Dimopoulos 1996) but no assump- 
tion about their order of magnitude is here adopted. 

Some important works have recently dealt with MHD 
in an expanding universe (Holcomb 1989, 1990; Dettmann, 
Frankel, & Kowalenko 1993; Gailis et al. 1994; Gailis, 
Frankel, & Dettmann 1995; Brandenburg, Enqvist, & Ole- 
sen 1996). However our objective is not MHD, but the in- 
fluence of magnetic fields on the formation of large scale 
structure. In these papers the metric is unperturbed. Here 
magnetic fields themselves are responsible for perturba- 
tions in the metric, which induce motions and density in- 
homogeneities, which in turn affect the perturbed metric, 
and possibly the magnetic fields. 



We have not included either protons or electrons in 
the system of equations. In a first attempt to solve this 
problem, this omission can be accepted, especially when 
we are considering an epoch of the Universe that is domi- 
nated by photons in which charges are considered to play 
a minor role. Nevertheless, W x B exists, which creates a 



charge current (see eq. (44)) and has an influence on the 
remaining equations. The existence of charges is implicitly 
assumed to support magnetic fields and electric currents, 
but equations of this third component are not explicitly 
considered. The infiuence of magnetic fields on the photon 
inhomogeneities lies in the fact that curvature is decided 
by the energy-momentum tensor, which is modified taking 
into account the magnetic contribution. This contribution 
is probably too small to affect the expansion itself, but 
not so as to have a large influence on the internal struc- 
ture. The contribution of magnetic fields to the energy- 
momentum tensor is very different with respect to the 
contribution of photon and baryon densities alone. This 
inclusion is not, therefore, trivial and indeed the results 
are clearly different. In the presence of magnetic fields, in- 
homogeneities evolve in a completely different way, as is 
demonstrated here in this particular epoch in the lifetime 
of the Universe. 

In agreement with the Cosmological Principle, we con- 
sider than no mean magnetic field exists at cosmological 
scales, so < B >= (or rather, we demonstrate in Ap- 
pendix (^) that < B >= in a Robertson- Walker met- 
ric) . Classically, this condition is equivalently reached from 
V • B = 0, Gauss theorem and the Cosmological Principle. 
However random magnetic fields do exist at lower scales in 
characteristic cells. These fluctuative magnetic fields, even 
with random orientations, are present everywhere, so that 
< > is non-vanishing. There is no mean magnetic field 
in the Universe, but there is a mean magnetic energy. No 
assumption about its value, such as that of equipartition 
with the CMBR energy density or any other hypothesis, 
is made "a priori" . 

In this paper we obtain the equations and derive the 
basic conclusions. In forthcoming papers we will deal with 
inhomogeneities affected by selected particular magnetic 
configurations, and with the infiuence of a large scale mag- 
netic field distribution on the large scale density distribu- 
tion. Paper II deals with the infiuence of magnetic fiux 
tubes on the distribution of the density, showing that pri- 
mordial magnetic fields can be responsible for the observed 
present filamentary structure 

The curvature has been set equal to zero, which is, in 
any case, a good assumption for this epoch. 

2. The mean magnetic field 

Before considering the interrelation between fluctuative 
magnetic fields and density inhomogeneities, let us con- 
sider the mean magnetic field and its influence on the mo- 
tion of the Universe as a whole. It is intuitive that an 
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isotropic universe cannot possess a mean magnetic field. 
Nevertheless, it is demonstrated in Appendix (^) that the 
existence of a mean magnetic field is incompatible with the 
Robertson- Walker metric, by examining the form adopted 
by Maxwell equations, the equation of motion and the Ein- 
stein field equations in a magnetized universe. Therefore, 
we consider < B >= 0. However, magnetic fields may be 
present (and actually they are) in smaller cells, with the 
direction of the field being random at larger scales, so that 
we assume < >^ 0. There is no mean magnetic field 
but there is a mean magnetic energy. 

In Einstein field equations we must include the mean 
magnetic energy-momentum tensor, deduced from: 
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where F^'' is the Faraday tensor. We adopt for the 
contravariant-covariant form of the Faraday tensor in a 
Minkowskian frame 
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where 



B = 



and E and B are the electric and magnetic three- vectors, 
as seen by an inertial observer. 

The Robertson- Walker metric for k = can be written 

as 



500 = -1 
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(5) 



We then obtain 
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Using the Robertson- Walker metric we obtain the other 
forms of the Faraday tensor 
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We then assume infinite conductivity, so that electri- 
cal fields in the rest frame of the charged particles van- 
ish. Magnetic fields are assumed to be tied to charges. 
Cheng & Olinto (1994) showed that the effects of finite 
conductivity in the early universe may be neglected. The 
electromagnetic momentum- energy tensor becomes 
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We also have 
< B >=0 



< Bi >-- 



<B'> 



--< Bl >=< Bl > 



(11) 



(12) 



(13) 



(14) 



(15) 



where (as usual) Latin indexes denote only spatial coor- 
dinates, and R is the cosmological scale factor. Through- 
out the paper we will consider that R is measured tak- 
ing its present value Rq as unity. Therefore, R is di- 
mensionless, i? — 1 at present, and we have approxi- 
mately R = z^^, with z being the redshift, taking into 
account that we are dealing with times long before Re- 
combination. Now we take into account the transforma- 
tion of coordinates that transforms — dia.g(— 1, 1, 1, 1) 
into Qfii, — diag{—l,R'^,R^,R^). This transformation is 
achieved by 



^ „ — -7^ ~ diaqfl, — , — , — ) 

" dx"" ' i?' i?' r' 



(6) 



The subindex M denotes "magnetic" 
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This is easily checked. For instance 47rT^} = R ^ < 
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and 
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is the magnetic energy density, we have as the equation of 
state 
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The form of the magnetic energy-momentum tensor is 
identical to the form of the radiative energy-momentum 
tensor. For photons we consider 
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where pr is the radiative hydrostatic pressure, cr the 
radiative energy density, and we consider = 1 and 
W = 0. Therefore 
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where 

e = Em + Efl 
and 

P = rM +Pr 



(22) 
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are the total density and hydrostatic pressure. Given the 
similarity of form between and T"^ the Einstein field 
equations must provide familiar results, i.e. the same ex- 
pansion and cooling laws that hold for a purely radiative 
universe. 

It is straightforwardly obtained for the Robertson- 
Walker metric that 
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with R = as usual, and therefore 



If Em <C e_R wc find the familiar result T cx R~^j but in 
general, we must specify a relation between cm and cr 
to deduce both. For instance, if there was at some epoch 
cquipartition of radiative and magnetic energy em = ^r, 
both would decrease at the same rate and the condition 
of cquipartition would be conserved. If this condition were 
confirmed at present (as suggested by Kronbcrg 1994) the 
epoch of the Universe dominated by photons would be- 
come co-dominated equally by photons and by magnetic 
fields. However, wc emphasize that this cquipartition con- 
dition is not assumed (or deduced) in this paper. If it held, 
the relation R oc would be maintained, but now the 
constant of proportionality would be higher, by a factor 
of less than y/2, and the expansion would be faster. 

If we assume that the magnetic energy density is much 
less than the radiative one, as is usually done, and indeed 
this is what our results suggest, then we would obtain the 
obvious result that the expansion is unaffected by mag- 
netic fields. 



3. The perturbed quantities 

As usual, any quantity is substituted by its mean value 
plus a fluctuating quantity, which in the linear approxi- 
mation is negligible, i.e. terms including products of two 
fluctuating quantities are considered second order terms. 
For instance €r ^ eR + Ser and pR pr + 6pR. As there 
is now no possibility of confusion we will eliminate the 
subindex R and write e + Se and p + Sp. 

By considering the transformation g^fj gap + 5gap 
we will call 6ga/3 = hap- But some transformations of 
the metric tensor do not mean real physical changes. 
As argued by Weinberg (1972), it is possible to choose 
hoi = hoo = 0. We benefit here from this choice. It is 
necessary to calculate /i*'^ defined as Sg"'^. The metric 
tensor must match (gajs + hap)(g^'' + h!^-"') = 5"^ hence 



+ haisg + gafjh^'^ = S2, neglecting higher order terms. 
Therefore h^'^ = ~g^^"hacrg"'' , and we have for each com- 

(25) ponent ht"^ = —R~^hi3^ again with ft,™ = = 0. 

hafj is equivalent to a three-dimension tensor, as any 
component containing the time subindex vanishes. When 

(26) using three-dimension formulae we term hij = h.lis trace 
hii plays an important role and will simply be called h. 
We also term ft* = h" = —R~'^hii = —R~^h. 

(27) W is the four-velocity of the photon fluid. We have 
5U'^ = 0, as f/" = 1 with no perturbation. When dealing 
with three-dimension formulae we term 6U^ = u. As we 
are using comoving coordinates, the unperturbed velocity 

(28) is U' = 0. 
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It is easily calculated that the components of the per- 
turbed affine connection SV^^^ vanish except 
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For quantities of electromagnetic nature, we would have 
for instance Bi — s- Bi + SBi. But as shown in the preceding 
paragraph the mean quantity is null and therefore we may 
use Bi instead of SBi . Therefore B , E 
perturbed quantities. 



Jem EM ^'"^ 



4. The perturbed Maxwell equations 

We will obtain the perturbed equations from Maxwell's, 
motion, energy and Einstein's field equations. Not all these 
equations are independent but they are all useful. The 
general procedure will be as usual. After substituting any 
generic quantity A by A + 5 A, the original equation is 
subtracted and second order terms neglected, because we 
only consider linear perturbations. As an exception, when 
perturbing Maxwell's equations we will not neglect second 
order terms, assuming that variations in the metric and in 
the Faraday tensor are uncorrelated. The reason is that, 
because the result is the same, and because some conclu- 
sions are very important in forthcoming paragraphs, we 
have preferred to be extremely cautious. 

Let us begin with the second set of Maxwell equations: 



this yields the following perturbed equation 
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Note that the six last terms on the left hand side are 
second order terms. 

In the absence of electric fields, from the [0,0,0], 
[0,0,?], [0,i,0] and [0, components, we obtain just 
= 0. From [0,i,j] 



df^'' ~ ^ 

and for [0, 1, 2] (for instance) 
d{R^B^) 



dt 
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From [i, j, k] 
_d_ _d_ 
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For [1,2,3] for example 

V • B = 

The whole set of equations is reduced to 

V • B = 



constant 



(35) 

(36) 

(37) 
(38) 



The first one is familiar. The second one is also famil- 
iar but not obvious. It brings to mind the condition of 
frozen-in magnetic field lines in present cosmic plasmas. 
What this equation shows is that the magnetic pattern 
always remains the same, and the magnetic strength is re- 
duced by just the effect of expansion. The original pattern 
is conserved, even if we are considering small motions of 
the photon fluid associated with density inhomogeneities. 
Neither do small metric fluctuations alter the magnetic 
pattern. This property is also obtained for the linear New- 
tonian post-Recombination epoch as shown in Appendix 

(i- 

From the first set of Maxwell equations 

J^"f„ = (39) 

with the above procedure and the same change of nomen- 
clature, we obtain 



_d_ 
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(40) 



The second and the fourth terms are second order terms, 
which were included in the derivation. One of them van- 
ishes, that is {6T^^)T°"' = 0. The resuh is 
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As [(V/i) X B] is negligible, we simply obtain that 
J° = (43) 
i.e. macroscopic neutrality, a familiar result, and that 
1 



J = 



47ri?2 



V X B 



(44) 



which simply provides J, the electric current, when B has 
been calculated. 
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5. The perturbed equation of motion-energy 

We start with the equation of motion 

r"f^ = 

where 
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which is the equation of energy conservation. For the ith 
component: 
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The electromagnetic contribution is as in equation ( |ll|) 
but now taking into account E = 0. We now obtain 

Ap^U" + + 
-pa ^cfi I p/5 ^aa , 
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However, it is possible to show that W ■ h — 0. This 
comes from the well known property of the tensor metric 
<?Q/3;7 = 0, and therefore gaf3;t3 — and gi/3-p = 0, which 
becomes in our case gip-p + hifs-f} where gi/3 now corre- 
sponds to the Robertson- Walker metric. We have 



dg 



'i/3 



" dxf^ 
hence 

hip-p = 

and 

U — - i pinal3 - i ppil-iij — 



dxi 



(51) 
(52) 
(53) 



With this simplification and some algebra, we have for the 
ith component 
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which is the equation of motion. 



6. The perturbed Einstein field equations 

We start with the Einstein field equations in the form 
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where T"^ is given by equation (46). For the electromag- 
netic contribution we have 
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R^i, is Ricci's tensor, where 
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as in Weinberg (1972). Then, we obtain 
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7. Density perturbations 

Let us define the relative perturbed density as: 

^ Sp Se 
P e 



(64) 



{6 is not so defined by other authors. Some others adopt 
6 ~ Se/{e + p) which becomes in this case 3/4 times the 
value of our 6). The final purpose is to deduce S{t,r), i.e. 
how density inhomogeneities evolve, under the action of 
gravity and an unperturbed underlying magnetic config- 
uration. To find the differential equation for S, we must 
perform some algebra because our set of equations is still 
rather large and complicated. We must first choose which 
equations will be handled, because as stated before, they 
are not all independent. Our unknowns are 6p, u, h, as 
B has already been determined as a function of time. 
We therefore have ten unknowns. The Einstein field equa- 
tions are also made up of ten equations, and from them 
the equations of motion and energy are derivable. On the 
other hand. Maxwell's equations have already been used 
to deduce the evolution of the magnetic field and the elec- 
trical current. Though we could therefore use only the 
Einstein field equations, the final objective is straightfor- 
wardly reached with the equation of motion (p4), the equa- 
tion of energy balance (49), and equations (|62|) and (63) 
from the whole set of Einstein field equations, selected be- 
cause they do not contain the tensor h but just its trace 
h. The unknowns are now 5p^ u and h. We still have more 
equations than unknowns. Let us rewrite the equations 
with the introduction of "present day" quantities, with 
the subindex 0. 



Po = pR 



Bq — BR^ 



(65) 



(66) 



ho = hR 



(68) 



Here po and B^ would actually be the present values 
of p and B if our equations were valid for all subsequent 
periods. The lawpi?^ — constant is valid for photons in the 
Acoustic and Post-Recombination eras. However B^ has 
recently been affected by non-linear effects and probably 
does not coincide with the present magnetic field. While po 
and Bo are constant, 5po and ho are not. Neither Spo nor 
ho represent actual present values, but they are preferred 
as variables because in the variations of Spo and /iq the 
effect of pure expansion is suppressed. We then have: 
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As suggested by the equations themselves, we now change 
the temporal and spatial variables, and use the definition 
( p4| ) of 5. As new spatial variables we will use 

x[ = ^x, = 2.73 X IQ-^^x^ (73) 

Re 

where K is the constant in the expansion law 

R^ = Kt (74) 

in the radiation-dominated era. Its value is 

is: = 4^/2^= 2.73 X 10"2°s"i (75) 

and Re is the cosmological scale factor at the last time 
in the period considered in this paper, which is close 
to Equality. For po we have adopted the value 8.84 x 
10-''2s-2_ should be noted that Xi are comovmg co- 
ordinates measured in seconds. They now coincide with 
present coordinates. (It is easily obtained that 1 Mpc = 
1.03 X lO^^'s). However x'^ is dimensionless 



x' = 0.28 



/ d 



\Mpc 



(76) 



where d is the distance in Mpc. 

As a temporal variable let us choose 



T = ln^ = 21n:^ 
t R 



(77) 
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and therefore 



8. Conclusions 



dt 



K 



(78) 



where te is the last time of the period considered here, 
close to Equality, corresponding to i? = 10~^, i.e. 3.66 x 
lO^s. r is a time variable increasing backwards from fu- 
ture to past. We then introduce some time independent 
functions defining the magnetic pattern: 



X 



247rpo 



Bo • V'Bo , OV7' V 
n = h o V A 



2m = 



-V' • -n-VX 
3 



(79) 



(80) 



(81) 



where V is the gradient using instead of Xi . In the unit 
system we are using 1 Gauss = 8.61 x 10~^^s~^. In this 
way all the magnitudes of the different quantities involved 
are close to unity and the equations become extremely 
simple. 

We also change the nomenclature, so that a point over 
a quantity's symbol now means its derivative with respect 
to r, instead of with respect to t. Thus we have 



2ReU - AReii + V'(5 + n = 



- 3(5 + 4i?ee"^V' • M - 2/io = 



2/io + 3^ + 3X = 



(82) 



(83) 



(84) 



(85) 



V'/l'o + 2ReU = 



Taking V' in equation ( ^4D and ^ in equation ( p5[ ) 
and subtracting, inserting this ii in equation (^) to ob- 
tain M, calculating ^ in equation (|8^), taking into account 
the obtained functions u and ii in equation (^), we ob- 
tain a differential equation containing only (5 as a variable, 
<5(t,xD 



2e- 







(86) 



This is our basic equation, an elliptic linear second order 
differential equation with variable coefficients. The com- 
plete discussion may have complications arising from its 
elliptic nature (if we prefer not to choose a boundary con- 
dition for T = 0) . The following section contains a prelim- 
inary insight. 



As previously mentioned, the implications of various indi- 
vidual magnetic patterns and of a cosmological magnetic 
spatial spectrum in the formation of structures will be 
studied in forthcoming papers. In particular, we will ex- 
amine in paper II the effects of a magnetic flux tube, the 
simplest structure with cylindrical symmetry. Our basic 
results will be confirmed with this simple particular case. 
We will now examine the general results obtainable from 
the differential equation ([8q). 

The primordial magnetic field spatial spectrum re- 
mains time independent in the linear approximation con- 
sidered here, being diluted just by the general expansion 
but conserving patterns and relative values. Therefore, 
present large scale magnetic field patterns in the inter- 
galactic medium, if we are eventually able to observe them, 
could reveal very early processes concerning magnetogene- 
sis. Complementary to the standard methods reviewed by 
Kronberg (1994), the measurement of the highest energy 
cosmic rays (Lee, Olinto, & Siegl, 1996) and other indi- 
rect methods (Plaga 1995; Kronberg 1995; Battaner et al. 
1991) are promising in the future. There is even the pos- 
sibility of measuring magnetic fields at the last scattering 
surface (Kosowsky & Loeb 1996). 

The fact that magnetic inhomogeneities are conserved 
until recent epochs in the history of the Universe supports 
the works by Wasserman(1978) and by Kim et al. (1994) 
in which a preexisting magnetic configuration determines 
the evolution of density inhomogeneities. The influence of 
this primordial magnetic field spectrum on the formation 
of large-scale structures, clusters and superclusters, as well 
as of galaxies, is very important. 

If magnetic fields have enough strength they can have 
a clear influence on the origin and evolution of density 
inhomogeneities. For very early times, with very large r, 
or equivalently for very large scale structures, the last two 
terms in eq. ( p6| ) would become smaller and the resulting 
equation would be easily integrated to give 



-X + cie + C2e7 



(87) 



This equation will be analyzed in detail in paper II. We 
now see that for very large r, b would become very large 
unless C2 = 0. For large r we would then have b = —X . 

Observe that when X — ^ this equation is just the 
solution in the absence of magnetic flelds, expressing the 
growth of the compressional mode as summarized in the 
classic book by Weinberg (1972), giving (5 cx t cx R? . 

Then, magnetic fields may provide an alternative 
mechanism to generate the promordial spectrum of den- 
sity inhomogeneities, which are then amplified by gravi- 
tational instabilities. Other primordial mechanisms, such 
as quantum fluctuations at Inflation, or exotic disconti- 
nuities produced by phase transitions, cannot be rejected, 
but magnetism is a very interesting additional possibility 
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as, if (5 were null initially, and only the homogeneous Uni- 
verse were perturbed by random magnetic fields, 6 would 
have been created. 

On the other hand, not only are primordial magnetic 
fields able to originate density structures at very early 
epochs, but they have a direct influence on the evolution 
of structures at R ~ 10^^i?o. The effect of this direct 
influence may be of a very different nature, depending on 
the initial conditions and on the type of the magnetic field 
pattern, producing concentrations of photons which would 
eventually be the site of baryon concentrations. 

The term X only depends on the magnetic field energy 
density, but the last term, including m, depends on the 
magnetic field as a vector, thus producing anisotropics in 
the evolution of density inhomogeneities. This anisotropic 
evolution would be more important at relative recent 
times, close to Equality, and for small scale structures. 
Later on, heat conduction and viscosity would eliminate 
small scale structures and, after Recombination, non- 
linear effects would complicate the initial structure, but 
the primordial distribution of magnetic fields and that of 
the density at the epoch of Equality would probably still 
be recognizable at present. 

It is important to obtain the order of magnitude of the 
field strength able to affect the evolution of the density in- 
homogeneities. The magnitude of X must be at least of 
the same order than S, which can have typical values of 
10^^ in the period of time considered. If X were unity, 
we would have equipartition between magnetic and radia- 
tive energy densities and we would have Bq « 3/iG. But 
if X w 10^^, we see from its definition {X depends on 
Bq) that Bq Ki lO^^G. Therefore, magnetism affects the 
evolution of the density pattern if Bq > lO^^G. Higher 
fields are not expected as they would produce a very 
fast density evolution, incompatible with present observa- 
tions. The magnetic field may not differ very much from 
Bo ~ IQ-^G. 

This is an equivalent-to-present value. Real magnetic 
fields at any time are calculated taking into account 
BR^ = Bq. It is important to calculate B at the epoch 
of nucleosynthesis. With R — 2.3 x 10~^ at the end of 
nucleosynthesis, we have Bn ~ 2 x lO^G, perfectly com- 
patible with the limits of Grasso and Rubinstein (1995) of 
Bat < 3 X 10i°G. 

A. Mean magnetic field in the Robertson- Walker 
metric 

The purpose here is to show that the mean magnetic field 
is null in a universe obeying the metric of Robertson- 
Walker, in order to justify this assumption in the paper. 
First consider the Maxwell equations: 



where now (in contrast with the nomenclature above used) 
J-"^ and Jf^ are the unperturbed Faraday tensor (a func- 
tion of the unperturbed mean electric and magnetic fields 
E and B) and the unperturbed charge density electric 
current vector. All these quantities have been assumed to 
be n ull throughout the paper. We would now have, from 







+ 



"yet 



T^cr •77 "per •77 "per •77 

From this equation we obtain 

— (i-i?)-2-£; 



V • £; = 

djR^Bs) d{R^E2) _ dEi 
dt dx^ dx^ 

(and other similar formulae) 
E = 

V • B = 



(A3) 

(A4) 
(A5) 

(A6) 

(A7) 
(A8) 



With equation ( |A7| ) in ( |A4| ), ( |A5| ) and ( |A6D we have 

BR"^ = constant (A9) 

Therefore the Maxwell equations under a Robertson- 
Walker metric are not incompatible with a mean magnetic 
field. The other set of Maxwell equations simply tells us 
that the Universe should be macroscopically neutral and 
that electric currents given by 



J = 



1 



V X S 



(AlO) 



exist. Neither is the law of conservation of momentum- 
energy very restrictive with respect to the existence of a 
cosmological magnetic field. We can similarly benefit from 
the derivation carried out in section (^ . Now, again, B is 
the unperturbed magnetic field. For energy conservation 
we would have 



de 
dt 



d_ 

dt 



B^ 



R Stt 



(All) 



Q/3 



-47rJ^ 



(Al) 



where e is the photon energy density. If a;^ is the direc- 
tion of the magnetic field, we have for the conservation of 
momentum, taking into account the homogeneity of the 
Universe -P^ = 



.B 



(A2) V— = V.(BB) 



(A12) 
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which does not imply B = at all. This result is however, 
deduced from the Einstein field equations. We again ben- 
efit from the above calculations, with B now being the 
unperturbed magnetic field. It is obtained that 





-{RR + 2R^)5z J 
'B^ 



= -2 




^R^BB 



' pRH, 



(A13) 



Components (1,2), (1,3) and (2,3) of this equation yields 



B1B2 = 



B1B3 ^ 



B2B3 = 



(A14) 



(A15) 



(A16) 



This does not mean that Bi = B2 = B3 = 0. One of 
them may be non- vanishing, for instance Bi. But then 
components (1, 1) and (2, 2) yield 



{RR + 2R^) = 2R^Bf - R^B{ - SirpR 



(A17) 



- {RR + 2R^) = -R^Bf - 87rpi?2 (a18) 
and subtracting we obtain Bi = 0, and therefore 



B = 



(A19) 



in a Robertson- Walker metric. 

We conclude that in a universe with a Robertson- 
Walker metric, the existence of B, even if compatible with 
the Maxwell equations and with the equations of conserva- 
tion of energy-momentum, is incompatible with the Ein- 
stein field equations. 

From (A13) it is then obtained 



R 

- + 8'Kp + B^ 
R 



r>2 

RR + 2i?2 = 87rpi?2 + —R^ 



(A20) 



(A21) 



Eliminating R in these two equations it is obtained 
87rei?3 1 



3 3 



^B'^R^ = R^R 



(A22) 



which is the generalization of the Einstein-Friedmann 
equation, and which was implicitly used when obtaining 

(il)- 



B. Linear Newtonian evolution of magnetic fields 

It will be shown that in the linear Post-Recombination 
epoch, the magnetic structures evolve according to the law 
3R^ — constant, with 3 being the amplitude of the mag- 
netic propagating inhomogeneity (we are not now using 
comoving coordinates) in a similar way to their evolution 
in the Radiation dominated epoch. 

Let us consider the induction equation 

dB 



— ^\/x(vxB) = B-\/v-v - VB - BV ■ v 

dt ^ ' 



(BI) 



As the mean field vanishes, let us again term 5B — *■ B. 
After the introduction of perturbations, this becomes 

dB 



dt 



B Vv-v ■ VB - BV • V 



(B2) 



where v is the mean velocity. We would have in general 
V ^ V + u, with u being the fluctuative velocity, but it 
would be present only in second order terms. Therefore v 
obeys 



R 

where r is the position vector. With 

V-u = -8 
R 

we obtain 

dB R „ R 

— = -2—B r VB 

dt R R 

The solution would be of the form 
B = BS 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



where 13 is the amplitude that only depends on time, and 
£ is a function only of r : 



(B8) 



where q is constant. Substituting in the induction equa- 
tion, we have 



R Q • V ■ O ' T • 

B -f 2-B - i^RB + i^RB = 
R rt^ R^ 

hence 

B + 2^B = 
R 

and 

BR^ = constant 



(B9) 



(BIO) 



(Bll) 



It can therefore be concluded that the magnetic spa- 
tial primordial pattern has been preserved until very re- 
cently when the back-reaction of the velocity field onto the 
magnetic field became important. The importance of the 
back-reaction has been pointed out by Kim et al. (1994). 
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